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Abstract 

A search for a Higgs boson in the four-lepton decay channel H — > ZZ, with each Z 
boson decaying to an electron or muon pair, is reported. The search covers Higgs 
boson mass hypotheses in the range 110 < mn < 600 GeV. The analysis uses data 
corresponding to an integrated luminosity of 4.7fb~ 1 recorded by the CMS detector in 
pp collisions at yjs = 7 TeV from the LHC Seventy-two events are observed with four- 
lepton invariant mass > 100 GeV (with thirteen below 160 GeV), while 67.1 ± 6.0 
(9.5 ± 1.3) events are expected from background. The four-lepton mass distribution 
is consistent with the expectation of standard model background production of ZZ 
pairs. Upper limits at 95% confidence level exclude the standard model Higgs boson 
in the ranges 134-158 GeV, 180-305 GeV, and 340^65 GeV. Small excesses of events 
are observed around masses of 119, 126, and 320 GeV, making the observed limits 
weaker than expected in the absence of a signal. 
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The standard model (SM) of electroweak interactions HHS relies on a scalar particle, the Higgs 
boson, associated with the field responsible for the spontaneous electroweak symmetry break- 
ing HH9|. The existence of the Higgs boson has yet to be established experimentally, while 
its mass, m^, is not fixed by the theory. Direct searches for the SM Higgs boson at the LEP 
e + e collider and the Tevatron pp collider have led, respectively, to a lower mass bound of 
m H > 114.4 GeV [10], and to an exclusion in the range 162-166 GeV El, at 95% CL. Indirect 
constraints from precision measurements favour the mass range mn < 158 GeV fT2\ at 95% CL. 
The inclusive Higgs boson production followed by the decay H — > ZZ is expected to be one 
of the main discovery channels at the CERN proton-proton (pp) Large Hadron Collider (LHC) 
for a wide range of wjh values. Using the H — > ZZ and the H — > WW decay channels, the AT- 
LAS collaboration have excluded at 95% CL the mass ranges 145-206 GeV, 214-224 GeV, and 
340^50 GeV mMa. 

In this Letter, an inclusive search in the four-lepton decay channel, H — > ZZ — > £ + £~ £' + £'~ with 
I, t! = e or }i, abbreviated as H — > 4£, is presented. The analysis is designed for a Higgs boson 
in the mass range 110 < win < 600 GeV. It uses pp data from the LHC collected at y/s = 7 TeV by 
the Compact Muon Solenoid (CMS) experiment during 2010 and 2011. The data corresponds 
to an integrated luminosity of 4.7 fb -1 . The search relies solely on the measurement of leptons, 
and the analysis achieves high lepton reconstruction, identification, and isolation efficiencies 
for a ZZ — > A£ system composed of two pairs of same-flavour and opposite-charge isolated 
leptons, e + e or in the measurement range > 100 GeV. One or both of the Z bosons 

can be off-shell. The background sources include an irreducible four-lepton contribution from 
direct ZZ production via qq annihilation and gg fusion. Reducible contributions arise from Zbb 
and tt where the final states contain two isolated leptons and two b jets producing secondary 
leptons. Additional background of instrumental nature arises from Z + jets events where jets 
are misidentified as leptons. 

Particles produced in the pp collisions are detected in the pseudorapidity range \rj\ < 5, where 
rj = — In tan(6V 2) and 6 is the polar angle with respect to the direction of the proton beam. The 
CMS detector comprises a superconducting solenoid, providing a uniform magnetic field of 
3.8 T in the bore, equipped with silicon pixel and strip tracking systems (\f]\ < 2.5) surrounded 
by a lead tungstate crystal electromagnetic calorimeter (ECAL) and a brass-scintillator hadronic 
calorimeter (HCAL) (\t]\ < 3.0). A steel /quartz-fiber Cherenkov calorimeter extends the cov- 
erage (\tj | < 5). The steel return yoke outside the solenoid is instrumented with gas detectors 
used to identify muons (\r/\ < 2.4). A detailed description of the detector is given in Ref. IITBl . 

Monte Carlo (MC) samples for the SM Higgs boson signal and for background processes are 
used to optimize the event selection and to evaluate the acceptance and systematic uncer- 
tainties. The Higgs boson signals from gluon-fusion (gg — > H), and vector-boson fusion 
(qq — > qqH), are generated with POWHEG [17J at next-to-leading-order (NLO) and a dedi- 
cated generator from Ref. ffl8]|. Additional samples of WH, ZH, and ttH events are generated 
with pythia [19] . Events at generator level are reweighted according to the total cross section 
c(pp — > H), which contains contributions from gluon fusion up to next- to-next-to-leading- 
order (NNLO) and next-next-to-leading-log taken from Refs. [20 -29] and from the weak-boson 
fusion contribution computed at NNLO in Refs. fl22l |3Q] - |34| . The total cross section is scaled 
by the branching fraction £>(H — > A£) ||22ll35ti38| . Interference effects in the 4e and 4fi chan- 
nels are calculated with PROPHECY4F (22H35H36|. The SM background contribution from ZZ 
production via qq is generated at NLO with POWHEG, while other di-boson processes (WW, 
WZ) are generated with pythia with cross sections rescaled to NLO predictions. The gg — > ZZ 
contribution is generated with GG2zz [39 | . The Zbb, Zee, Z7, and Z + light jets samples are 
generated with MadGraph [40] with cross sections rescaled to NNLO prediction for inclusive 
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Z production. The tt events are generated at NLO with POWHEG. The generation takes into 
account the internal initial state and final state radiation (FSR) effects which can lead to the 
presence of additional hard photons in an event. All events are processed through a detailed 
simulation of the CMS detector based on GEANT4 |4T| and are reconstructed with the same 
algorithms that are used for data. 

Collision events are selected by the trigger system that requires the presence of a pair of elec- 
trons (a pair of muons) with transverse energy (transverse momenta) for the first and second 
lepton above 17 and 8 GeV respectively. The trigger efficiency within the acceptance of this 
analysis is greater than 99% for signal in the 4e and 4p channels, and rises from about 97.5% 
at 7ft H = 120 GeV to above 99% at m H > 140 GeV in the 2e2p channel, within the acceptance of 
this analysis. 

Electrons are reconstructed within the geometrical acceptance, \ t] e \ < 2.5, and with p| > 7 GeV, 
by combining information from the ECAL and inner tracker l|42l |43 | . Electron identification 
selection requirements rely on electromagnetic shower-shape observables and on observables 
combining tracker and calorimetry information. The selection criteria depend on p| and \rj e \, 
and on a categorization according to observables sensitive to the amount of bremsstrahlung 
emitted along the trajectory in the inner tracker. Muons are reconstructed [44] within \ t]P \ < 2.4 
and Pj > 5 GeV, using information from both the inner tracker and the muon spectrometer. The 
inner track is required to be composed of more than 10 tracker-layer hits to ensure a precise 
measurement of the momentum. The efficiencies are measured in data, using a tag-and-probe 
technique [45] based on an inclusive sample of Z events. The measurements are performed in 
several ranges in p j and 1 t] \ . The product of reconstruction and identification efficiencies for 
electrons in the ECAL barrel (endcaps) vary from about 68% (62%) for 7 < p| < 10 GeV to 
82% (74%) at pf ~ 10 GeV, and reach 90% (89%) for p\ ~ 20 GeV. It drops to about 85% in 
the transition region, 1.44 < \r\\ < 1.57, between the ECAL barrel and endcaps. The muons 
are reconstructed and identified with efficiencies above ~98%. Lepton candidates are defined 
with a loose constraint on their isolation, by requiring the sum of the transverse momenta of 
tracks i within a cone around the lepton £ of AR = \J {r\ l — rj 1 ) 1 + (<p^ — (p 1 ) 2 < 0.3, where (p is 
the azimuthal angle, to have Ya Pj track/ Pr < 0.7. The lepton isolation efficiency for identified 
leptons with this very loose definition of isolation is found to be greater than 99%. 

We first require a Z candidate formed with a pair of lepton candidates satisfying 50 < m\ f % < 
120 GeV, pj > 20 GeV, and p T 2 > 10 GeV. The px thresholds ensure that the leptons are on the 
high-efficiency plateau for the trigger. The lepton pair is required to be well isolated using a 
combination of the tracker, ECAL and HCAL information. The sum of the combined relative 
isolation R[ S0 for the two leptons is required to satisfy R^ + R[ s 2 < 0.35, where for each lepton, 

Rise = (l/p T ) x (j^i P T/ track + £ t,ecal + 5jfc £ t,hcal)' with sums running over the charged 
tracks i, and the Ej from energy deposits in cells j and k of the ECAL and HCAL within a cone 
of radius AR < 0.3, respectively. The footprint of the lepton is removed from the isolation sum. 
The combined isolation efficiencies measured with data using the tag-and-probe technique are 
found to be > 99% for muons and between 94% to 99% for electrons. The isolation is made 
largely insensitive to the number of overlapping pp interactions by correcting for the average 
energy flow [46] per unit area measured as a function of the number of primary vertices. The 
ratio of the efficiencies measured with data and with simulated Z — > ££ events is found to be 
consistent with unity. The significance of the signed impact parameter of each lepton relative 
to the event vertex, SIP3D = where IP is the impact parameter in three dimensions and cr IP 
the associated uncertainty, is required to satisfy |SIP3d| < 4. The £ + £~ pair with reconstructed 
mass closest to the nominal Z boson mass is retained and denoted Z%. The Zi+X dataset thus 
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defined contains the samples used to estimate reducible and instrumental backgrounds, as 
well as the ZZ rates. In the next step, a subset of events is identified with at least a third 
lepton candidate. The Zi + £ events are used to measure misidentified lepton rates. A subset 
of events with at least a fourth lepton candidate of any flavour or charge is then identified. 
Together, the Zj + £ and Zj + ££' samples provide ways to estimate the remaining reducible 
(Zbb, tt) and instrumental (Z + light jets) backgrounds. For the signal, we select a second lepton 
pair, denoted Z2, from the remaining same flavour £ + £~ combinations, by requiring mz 2 > 
12 GeV, with the restriction > 100 GeV. For the 4e and Afi final states, at least three of 
the four combinations of opposite-sign pairs must satisfy tnu > 12 GeV. If more than one Z2 
candidate satisfies all criteria, the ambiguity is resolved by choosing the leptons of highest pj. 
The isolation and impact parameter are used to further suppress the remaining backgrounds. 
We require for any combination of two leptons i and /, irrespective of flavour or charge, that 
^iscy + Ki S0/ y < 0.35 and also impose | SBP30 1 < 4 for each of the four leptons. 

Finally, to select the four-lepton signal candidates, we require that the Z\ and Z2 masses satisfy 
m™" < m Zj < 120 GeV and m™ in < m Z2 < 120 GeV, with (m™ n ,mf^) = (50, 12) GeV defining 

the baseline selection and (m^™, mg™) = (60,60) GeV defining the high-mass selection. The 
baseline selection is used to search for the Higgs boson, and the high-mass selection is used to 
measure the ZZ cross section. 

The event yields are found to be in good agreement with the MC background expectation 
at each step of event selection. The ZZ and Z+X backgrounds dominate after the full event 
selection. The overall signal detection efficiency for the 4e (4fi, 2e2]d) channel is evaluated by 
MC and increases from «21% (59%, 35%) at m H = 120 GeV to ^35% (71%, 50%) at m H = 
140 GeV, reaching a plateau at «51% (81%, 63%) at m H = 200 GeV, and then slowly rising to 
~60% (83%, 72%) at mn = 350 GeV. The relative mass resolution estimated from MC signal 
samples is about 2.1% (1.1%, 1.6%) for 4e (4f/, 2e2^). 

The small number of observed events precludes a precise direct evaluation of background by 
extrapolating from mass sidebands. Instead, we rely on MC to evaluate the number of events 
expected from the ZZ background. The cross section for ZZ production at NLO, through the 
dominant process of qq annihilation and through gg fusion, is calculated with MCFM |47l - l49l . 
The theoretical uncertainties are computed as a function of m^, varying both the QCD renor- 
malization and factorization scales and the parton distribution functions (PDF) set following 
the PDF4LHC recommendations [50-54]. The uncertainties for the QCD and PDF scales for 
each final state are on average 8%. The number of predicted ZZ — > A£ events and their uncer- 
tainties after the baseline selection are given in Table [I] As a consistency check, an evaluation 
is made based on a normalization to the measured inclusive single-Z production, a procedure 
discussed in Refs. [55, 56 1. The measured rate of single Z bosons defined in this analysis is used 
to predict the total ZZ rate; making use of the ratio of the theoretical cross sections for ZZ and 
Z production, and the ratio of the reconstruction and selection efficiencies for the four-lepton 
and two-lepton final states. The results are in agreement with the ZZ rates reported in Table [I] 
within uncertainties. 

To estimate the reducible (Zbb, tt) and instrumental (Z + light jets) backgrounds, a region well 
separated from the signal region is defined by relaxing and inverting some selection criteria and 
verifying that the event rates change according to MC expectation. The event rates measured in 
the background control region are then extrapolated to the signal region. The control region for 
Z+X, where X stands for bb, cc, gluon or light quark jets, is obtained by relaxing the isolation 
and identification criteria for two additional reconstructed lepton objects (a measured track for 
muons, or a combination of a track and a cluster of ECAL energy deposits for electrons) indi- 
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cated as ^reco^reco- The additional pair of leptons must have like sign charge (to avoid signal 
contamination) and same flavour (e ± e ± , }i ± pi ± ), a reconstructed invariant mass mz 2 either satis- 
fying the baseline selection or the high-mass selection, and > 100 GeV. A sample Z\ + 1 xeco , 
with at least one reconstructed lepton object, is also defined for the measurement of the lepton 
misidentification probability, the probability for a reconstructed object to pass the isolation and 
identification requirements. The contamination from WZ in this set of events is suppressed by 
requiring that the imbalance of the measured energy deposition in the transverse plane is be- 
low 25 GeV. From the Z + 4eco^reco sample the expected number of Z+X background events in 
the signal region is obtained by taking into account the lepton misidentification probability for 
each of the two additional leptons. The number of background events expected in the signal 
region, normalized to the integrated luminosity, and the associated systematic uncertainties, 
are given in Table [ljfor the baseline selection in the range 100 < < 600 GeV. The reducible 
and instrumental background is found to be dominated by Z + light jets. A small residual 
contamination of Zbb remains at low mass while for the high-mass selection these reducible 
backgrounds are an order of magnitude smaller and therefore can be neglected. This was ver- 
ified by performing a measurement of Zbb and tt rates in a dedicated four-lepton background 
control region, defined by requiring a Z\ and two additional leptons satisfying an inverted 
SIP3D requirement, namely |SIP3d| > 5, and with relaxed isolation, charge, and flavour re- 
quirements. This ensures a negligible Z + light jets contribution in the four-lepton background 
control region, while the signal and the ZZ background are absent. To extract background rates, 
the reconstructed Z\ mass for the sum of the Zi + 2e, Z\ + 2}i, and Z\ + e]i final states is fit with 
a Breit-Wigner function convoluted with a Crystal Ball function [57| for the Z\ peak from Zbb 
and Chebychev polynomials for the description of the tt continuum. The extrapolation to the 
signal region relies on knowledge of, and the distinct features of, the SIP3D distributions for the 
Z2 leptons of the tt and Zbb backgrounds. The result is found to be compatible with the MC 
expectation in the signal region within the systematic uncertainty of 20%. 

Systematic uncertainties are evaluated from data for trigger efficiency (1.5%), lepton recon- 
struction and identification (2-3%), and isolation efficiencies (2%). Systematic uncertainties 
on energy-momentum calibration (0.5%), and energy resolution are accounted for by their ef- 
fect on the reconstructed mass distributions. The effect of the energy resolution uncertainties 
is taken into account by introducing a 30% uncertainty on the width of the signal mass peak. 
Additional systematic uncertainties arise from limited statistics in the reducible background 
control regions. All reducible and instrumental background sources are derived from control 
regions, and the comparison of data with the background expectation in the signal region is 
independent of the uncertainty on the LHC integrated luminosity of the data sample. This 
uncertainty (4.5%) 11581 enters the evaluation of the ZZ background and in the calculation of 
the cross section limit through the normalization of the signal. Systematic uncertainties on the 
Higgs boson cross section (17-20%) and branching fraction (2%) are taken from Ref. |22| . 

Recent studies f22l l59l l60l show that current Monte Carlo simulations do not describe the ex- 
pected Higgs boson mass line shape above ~300 GeV. These effects are estimated to amount to 
an additional uncertainty on the theoretical cross section, and hence on the limits, of about 4% 
at m H = 300 GeV and 10-30% for m H of 400-600 GeV. 

The number of candidates observed, as well as the estimated background in the signal region, 
are reported in Table [ljfor the baseline selection. The reconstructed four-lepton invariant mass 
distribution for the combined 4e, and 2e2}i channels with the baseline selection is shown in 
Fig.jl^ and compared to expectations from the backgrounds. The shape of the mass distribution 
below ffln = 180 GeV reflects the shape of the dominant qq annihilation process [61 1. The low 
mass range is shown in Fig. lib together with the mass of each candidate and its uncertainty. 
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Figure 1: a) Distribution of the four-lepton reconstructed mass for the sum of the 4e, 4^, and 
2e2}i channels, b) Expansion of the low mass range with existing exclusion limits at 95% CL; 
also shown are the central values and individual candidate mass measurement uncertainties. 
Points represent the data, shaded histograms represent the background and unshaded his- 
togram the signal expectations. 

The reducible and instrumental background rates are small. These rates have been obtained 
from data and the corresponding distributions are obtained from MC samples. 

The measured distribution is compatible with the expectation from SM direct production of 
ZZ pairs. We observe 72 candidates, 12 in 4e, 23 in and 37 in 262^, while 67.1 ± 6.0 events 
are expected from standard model background processes. No hard photon (pj > 5 GeV) was 
found, outside the isolation veto cone that surrounds each lepton, that could be unambiguously 
identified as FSR. Thirteen candidates are observed within 100 < < 160 GeV while 9.5 ±1.3 
background events are expected. We observe 53 candidates for the high-mass selection com- 
pared to an expectation of 51.3 ± 4.6 events from background. This high-mass event selection 
is used to provide a measurement of the total cross section c(pp — y ZZ + X) x B(ZZ — > 4£) = 
28.1^4Q(stat.) ± 1.2(syst.) ± 1.3(lumi.) fb. The measurement agrees with the SM prediction at 
NLO [47J of 27.9 ± 1.9 fb. The local p-values, representing the significance of local excesses 
relative to the standard model expectation, are shown as a function of mn in Fig. [2^, obtained 
either taking into account or not the individual candidate mass measurement uncertainties, 

Table 1: The number of candidates observed, compared to background and signal rates for 
each final state for 100 < tn^i < 600 GeV for the baseline selection. For the Z+X background, 
the estimations are based on data 



Channel 


4e 


A}i 


lelji 


ZZ background 


12.27 ± 1.16 


19.11 ± 1.75 


30.25 ± 2.78 


Z+X 


1.67 ± 0.55 


1.13 ± 0.55 


2.71 ± 0.96 


All background 


13.94 ± 1.28 


20.24 ± 1.83 


32.96 ± 2.94 


m H = 120 GeV 


0.25 


0.62 


0.68 


m H = 140 GeV 


1.32 


2.48 


3.37 


m H = 350 GeV 


1.95 


2.61 


4.64 


Observed 


12 


23 


37 
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for the combination of the three channels. Excesses are observed for masses near 119 GeV and 
320 GeV. The small ~2a excess near 320 GeV includes three events with > 50 GeV. The 
most significant excess near 119 GeV corresponds to about 2.5c significance. The significance 
is less than 1.0c (about 1.6c) when the look-elsewhere effect [62 1 is accounted for over the full 
mass range (for the low-mass range 100 < < 160 GeV). The local significances change 
only slightly when including candidate mass uncertainties, instead of using the average mass 
resolution, e.g. rising to 2.7a around 119 GeV and reaching 1.5c around 126 GeV. 




1.4 1 I i i I i i i I I I i i i i I i i i i I I 1 1 1 — i — 

110 200 300 400 500 600 110 200 300 400 500 600 



m H [GeV] m H [GeV] 

Figure 2: a) The significance of the local excesses with respect to the standard model expectation 
as a function of the Higgs boson mass, without (blue) or with (red) individual candidate mass 
measurement uncertainties, b) The observed and the median expected upper limits at 95% CL 
on c(pp — > H + X) x £>(ZZ — > Ai), normalized to the standard model cross section values a$M, 
for a Higgs boson in the mass range 110-600 GeV, using the CL S approach. The insets expand 
the low mass range. 



In absence of a significant clustering of candidates at any given mass, we derive exclusion lim- 
its. The exclusion limits for a SM-like Higgs boson are computed for a large number of mass 
points in the mass range 110-600 GeV, using the predicted signal and background mass distri- 
bution shapes. The choice of the step size in the scan between Higgs mass hypotheses is driven 
by either detector resolution, or the natural width of the Higgs boson. The signal mass distri- 
butions shapes are determined using simulated samples for 27 values of covering the full 
mass range. The shapes are fit using a function obtained from a convolution of a Breit-Wigner 
probability density function to describe the theoretical resonance line shape and a Crystal Ball 
function to account for the detector effects. The parameters of the Crystal Ball function are 
interpolated for the ?«h points where there is no simulated sample available. The shapes of 
the background mass distributions are determined by fits to the simulated sample of events, 
while the normalization is taken from estimates of overall event yields as described above. For 
each mass hypothesis, we perform an unbinned likelihood fit using the statistical approach 
discussed in Ref. [63 J. We account for systematic uncertainties in the form of nuisance parame- 
ters with a log-normal probability density function. The observed and median expected upper 
limits on c(pp -> H + X) x B{H -» ZZ) x B{ZZ -» U) at 95% CL are shown in Fig. |2j?. The 
limits are calculated relative to the expected SM Higgs boson cross section values asm, using 
the modified frequentist method CL S AMI l65ll . The bands represent the la and 2a probabil- 
ity intervals around the expected limit. These upper limits exclude the standard model Higgs 
boson at 95% CL in the m H ranges 134-158 GeV, 180-305 GeV and 340-465 GeV. The limits re- 
flect the dependence of the branching ratio £>(H — > ZZ) on niu- The worsening of the limits 
at high mass arises from the decreasing cross section for the H — > A£ signal. By virtue of the 
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excellent mass resolution and low background, the structure in the measured limits follows the 
fluctuations of the number of observed events. 

In summary, a search for the standard model Higgs boson has been presented in the four- 
lepton decay modes. Upper limits at 95% confidence level exclude the Higgs boson mass ranges 
134-158 GeV, 180-305 GeV, and 340^65 GeV. A major fraction of the explored mass range is 
thus excluded at 95% CL and the exclusion limits extend beyond the sensitivity of previous 
collider experiments. Excesses of events are observed at the low end of the explored mass 
range, around masses of 119 and 126 GeV, and at high mass around 320 GeV. These excesses, 
although not statistically significant, make the observed limits weaker than expected in the 
absence of a signal. At low mass, only the region 114.4 < nin < 134 GeV remains consistent 
with the expectation for the standard model Higgs boson production. 
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